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We investigate the pure penguin process B~ — > ir~ <f> using QCD factorization approach to calculate 
hadronic matrix elements to the a 3 order in some well-known NP models. It is shown that the 
NP contributions in R-parity conserved SUSY models and 2HDMs are not enough to saturate the 
experimental upper bounds for B — > <jyn. We have shown that the flavor changing Z 1 models 
can make the branching ratios of B — > <j)n to saturate the bound under all relevant experimental 
constraints. 



PACS numbers: 13.25.Hw, 12.38.Bx 

The process B — > <f>n, one of charmless two-body nonlcptonic decays of B mesons, is interesting 
because it is a pure penguin process. In particular, there are no annihilation diagram contributions for 
which results obtained with different methods are quite different [l[ 0] ■ Therefore the calculations of 
the hadronic matrix elements relevant to the process are relatively reliable because of no contributions 
coming from diagrams of annihilation topology. It proceeds through b — ► dss at the quark level, which 
is a b — + d flavor changing neutral current (FCNC) process. It is sensitive to new physics (NP) because 
all contributions arise from the penguin diagrams in the standard model (SM). 

The BaBar collaboration has recently reported the results of search for B — ► 4>tt 

Br(B° — > 07T°) = (0.12 ± 0.13) x 10" 6 , (1) 

< 0.28 x 10~ 6 , 
Br(B+ -> 07T+) < 0.24 x 10~ 6 , 

which enhance the precision of measurements but still are of the same order of mag nitude, i.e., O(10" 7 ), 
comparing with the previous results [4| 

Br(B° -» 0tt°) = (0.2±g;| ± 0.1) x lO" 6 , (2) 

< (1.2 ±0.8) x 10~ 6 , 
Br(B + -> </.tt + ) < 0.41 x 10~ 6 . 



The SM predictions for these decays modes in both BBNS (QCDF) and Li ct al. (PQCD) methods 
have been given (sj and the results given in the literature are much smaller than the data although there 
are significant disagreements among the literature. The twist-3 contributions to the decays are predicted 
to be small and Br(B + — > <f>ir + ) = (2 — 6) x 10~ 9 is given by using the QCDF method improved in 
calculating integrals which contain endpoint singularities, including the twist-3 contributions [5j|. The 
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agreed results are obtained in QCDF without including the twist-3 contributions in ref. [6j. The theo- 
retical uncertainty in treating integrals which contain endpoint singularities is roughly 30%. Therefore, 
there is a room for new physics (NP). 

The theoretical study for these decay modes in models beyond the standard model has also been 
done in a number of papers 0, @, 0|- The branching ratio for Br(B + — > (f>n + ) has been calculated in 
constrained minimal supersymmetric standard model (CMSSM) without imposing the constraint from 
B s — > 5]. In ref.p| calculations are done in MSSM, the topcolor assisted technicolor model 

(TC2), and the model with an extra vector like down quark (VLDQ) respectively. Only for VLDQ 
Br(B° — > 4>ir°) ~ 10~ 7 can be obtained. The analysis in MSSM uses the values of the mass insertion 
parameters which are taken from the paper in 1996 [8| and does not consider the contributions from 
neutral Higgs boson induced operators. In ref. [7] the upper bound of Br(B + — ► 07r + ) in 2002 Q is used 
to constrain parameters in R-parity violating supersymmetric models. 

In the letter we show that O(10~ 7 ) branching ratio for Br(B — > <\m) can not be reached in MSSM 
and two Higgs doublet models (2HDM) I, II, and III when all relevant constraints from experiments are 
imposed. It can be reached in a flavor changing Z' model under all relevant experimental constriants. 

To begin with, we record the SM analysis in ref.Q briefly. The AB = 1 effective weak Hamiltonian 
in SM is given by 



njSM 



G 



C 2 Q 2 



=3. ...,10 



Ci Qi + CV 7 Q77 + Cgg Q 



h.c. 



(3) 



where X p = V p bV* d , Q\ 2 ar e the left-handed current-current operators arising from IV-boson exchange, 
Q3.....6 and Q7,...,io are QCD and electroweak penguin operators, and and Q$ g are the electromag- 
netic and chromomagnetic dipole operators, respectively. Their explicit expressions can be found in, 
e.g., Ref. 10j. Due to the flavor and color structures of the final state (fm, the chromomagnetic dipole 



operator Qs g does not contribute to the decays, and the tree operators Qi : 2 contribute only through 
electromagnetic corrections which is small numerically. In QCDF up to the order of a s , the decay 
amplitude for B~ — > ir~(f> is [5| 



A(B~ -> TT-(j)) 

V2A(B° -> ttV) 



(4) 



where 



a 5 (^) = C 5 + lc 6 + ^^C 6 (-^-12), 



a?(?r0) 
af(7T0) 



C7 + § + S^ 1 Cs { - F - 12) + + 3C ^ 

C9 + 4<?io + ^^C W F+ ?PPUCi + 3C 2 ) , 

97T 



(5) 
(6) 
(7) 
(8) 



N 4vr N 

Due to the almost completely cancellations of the two terms in the Wilson coefficient combinations, 
C3 + C^/Nc and C5 + Cq/N c , ag is the biggest among a,'s, i.e., the contributions of the electroweak 
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penguin operators dominate. With values of input parameters in ref.jjjj], Br(B + — > 07r + ) = 4.5 x 10 
is obtained. 



We now turn to MSSM and 2HDM. The effective Hamiltonian in 2HDM and MSSM can be written 



as 



B0 



njnew 



niSM 

Gp 
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n/new 



Eh'4 E [^ft+ci^] 



i=ll,...,16 



(9) 



(10) 



where Qn to Qi6 are the neutral Higgs penguin operators and their explicit forms can be found in 
Rcf. 0, 11 1 with the substitution s — > d. The primed operators, the counterpart of the unprimed 
operators, are obtained by replacing the chirality in the corresponding unprimed operators with opposite 
ones. In ref.0 the neutral Higgs penguin operators Q,, i = 11, 16 are not considered. 



The contributions of new operators to the decay amplitude for B 

A new (B- -^TT-(j)) 
= V2A new (B° ^tt» 



IS 



F»-*« (mi) 2e*- PB , 



(11) 



where we have neglected 0(— £) terms, 



m B J 4, 
4e- PB f<f, 



(12) 



a'^i = 3,5,7,9, come from the contributions of the primed operators and their expressions can be 
obtained by replacing the Wilson coefficients in the corresponding unprimed operators with primed 
ones, <zi;li3 come from the contributions of neutral Higgs boson induced operators which we shall 
discuss later on. 

In MSSM new contributions to a^i = 3,5 and the contributions to a[,i = 3,5 come from the 
gluino-sbottom loop and depend on the mass insertion parameters (<^4)i3)^4 = L,R. The constraints 
on (Sa A )i3, A — L,R from the mass difference AMd have been reanalyzed, using NLO QCD corrections 
of Wilson coefficients and recent lattice calculations of the B parameters for hadronic matrix elements 
of local operators [l^] • Depending on the average mass of squarks and the gluino mass as well the CKM 
angle 7, (^,4)13, A — L,R, is constrained to be O(10 -2 ) in the (<5^i)i3 = (,5rr)i3 case which is assumed 
in ref. [6j. So in this case the new contributions to a^, i = 3, 5 and the contributions to a' i: i — 3, 5 are 
negligibly small, compared with the SM. Only for the case of single (5^ £ )i3 (or (6^)13) non zero, 
(^Li)i3 ( or {&Rn)ia) can reacn O(10 _1 ) under the constraint from the mass difference AMd [12j and 
consequently the SUSY contributions can make the branching ratios of B — > 4>ir reached O(10 -8 ). 



The contributions of neutral Higgs boson induced operators, 1111,13, & r e as follows. 



Oil, 13 



"4^ ~N~ 



(fs+f")C Q 



(13) 
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where 



du [ In 2 u + 2 In u — 21A 2 (u)] S (u) 



du 



dk / <i& In 



1 4k 2 



>"T. 



+ u 2 + u 



J (bk)V s (u,b) 



(14) 



, j _ 2ir tub U/b 

Is 



pB^i 



"I . 



\du] 



6 2 (b 1 + b 2 )bP B (^b)V s (v,b 2 ) 



x [fi p (u + v) V p (u, h) +m B (£,- v) V{u, h)] 



6(b 2 - b)I (by/u£mB) K (b^u^mB) + 0(b - b 2 ) IoW«("Ib) K (y / -uvmbb 2 ) 



(15) 



Note that there are no contributions of neutral Higgs boson induced operators at the leading order in 
the a s expansion in the approximation omitting 0{m 2 ,/m 2 B ) terms because of their Dirac structure. 

In MSSM and 2HDM model I, II, c['l 13 (m w ) ~ 0.037 because of the constraint from B d -> 
[13]. So the a s corrections from neutral Higgs boson induced operators are negligible. 

In model III 2HDM there are neutral Higgs boson-mediated FCNC at the tree level 14, 15, 
The Yukawa Lagrangian for quarks can be written as [17j 



-UM V U - DM D D 



2M W 



(H u cos a - h" sin a) UM V U + DM D D 



19 G° (UMm 6 U - DM Dl b D) 



2M W 



-JL—G+UVckm 
V2M W 

H° sin a + h° cos a 



Mu^{1+^)-Md\{1-1 5 ) 
M D i(l + 7 5 )-M t7 i(l- 7 5 ) 



U 
D 



V2 



D 



+ ~2 

- H+U 

- H~D 



1 



u ( e^i + 7 b ) - e^a - 7 & ) ] u - d[ + 7 b ) - e" T ~(i - 7 



1 



D 



V C k»lZ d \(1 + 7 5 ) ~ ^ckm~(1 - 7 5 ) 



1 



1 



D 



(16) 



where U and D represents the mass eigenstates of u,c,t quarks and d,s,b quarks respectively, the 
matrices t; D ' U are in general non-diagonal which parameterize the couplings of Higgs to quarks. The 
Yukawa Lagrangian for leptons can similarly be written. 

With the above Lagrangian, the Wilson coefficients of Qi,i = 11, 16, relevant to b — > dss are 
easily obtained: 



Cli,i3(myr) = 



sin 2 a cos 2 a 



G F (-X t ) y 8 



1 sTO b 



■5 



Cl2,14,15,16(mw) — 0, 

Cll,13 = Cll,13, Cl2, 14,15, 16( TO w) = > ( 17 ) 

where we have assumed £y is symmetric and real for simplicity. The parameters £bd,ss of the model 
can be estimated from experiments. We can extract £bd from the mass difference AMB d of neutral Bd 
mesons. In order to determine £ ss we use the data of r — > fifi + (JT , -B<j — > /i + /-t~ and r — > ^.P{P = 
7T°, T),rf). The parameter 6,^ has been estimated from the measured £?<2 — mass difference, and 
£dt < 7.3 x 10 -6 is give n [181 ] . We reanalyze the mass difference and get £db < 0.8 x 10~ 4 which is much 
larger than that in [18|. With the latest limits from experiments, Br(Bd — ► /i/i) < 3.9 x 1CP 8 [19( and 



_Br(r — * fififi) < 1.9 x 10 7 20| at 90% C.L., we obtain the updated results of the upper bounds of £ 



and („ T , {u U < 0.34 and < 0.004. Finally we can determine the bound on £ ss from r — > /j,tt°(t], rf) 



decays [2l| and the result is £ ss < 0.26. Using the bounds of ^6 an d £ ss , we deduce that Ci 3 (m w ) can 
roughly reach 0.1 with wia = mn = 300GeV, mh — !20GeV, which leads to that the a s corrections 
from neutral Higgs boson induced operators can contribute to the branching ratio of Br(B + — > 07r + ) 
by about 2.0 x 10~ 10 which is much smaller than that in SM and consequently negligible. 

Finally we consider flavor changing Z' Models. Because the contributions of electroweak operators 
dominate in SM it is expected that the new contributions from flavor changing Z' models in which Z' 
mediates vector and axil vector interactions would enhance the branching ratios of B — > fin significantly. 
The flavor changing Z' Models have been extensively studied f^ . For our purpose we write the Z' 
interaction Lagrangian in the gauge basis as 

C z ' = -g'^Z"* , (18) 
J p = X)^^ [( e V>x-)y p i + (Hn)i] P R\ > ( 19 ) 

id 



where g' is the gauge coupling constant of the U(l)' group at the Mw scale 1 , the sum extends over all 
fermion fields in the SM. There are in general FCNCs at the tree level in the Lagrangian (jT8j) . In the 
mass eigenstate basis, the chiral Z' coupling matrices are respectively given by 

= V ux e ux V^ x , B* = V dx e dx V} x , (X = L,R) (20) 

where V ux and V dx are the transformation matrices which make up-type quarks and down-type quarks 
into the mass eigenstates from the gauge eigenstates respectively. The usual CKM matrix is given by 
Vckm — V Ul VJ l . We assume that the first two generation diagonal elements of B^(q — u, ,d) are 



equal: B^ d = B* [24l27|. The effective Hamiltonian of the b — > dqq transitions due to the Z' mediation 



is 



H% = ^(^^) 2 i?^(M) y _ A ^(i?^(^) y _ A + i3«(^) y+A )+h.c. , (21) 



y/2 \gzM z , 

where Mz' is the mass of Z' which is larger than 850(76 V for g f ~ gz [2^ . 



sin tiw cos t>w 



1 We neglect the renormalization running effects from the ML scale to the Myy scale through the paper due to uncertainties 
in the parameters of the models. 
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Thus, the Z' mediated FCNC interaction would induce the color singlet QCD and electroweak penguin 
operators. It is straightforward to get the contributions to the Wilson coefficients of those operators 
from above effective Hamiltonian, (|21|) . Then the effective Hamiltonian relevant to B — > 4>ir in a Z' 
model can be written as 



G, 



(23) 



i=3,5,7,9 



where A t = VIVu, and 



zi 

3(5) 



3A 4 \ gi M z , 



z , 4 / g 2 M Z \ nL* / r>L(R) r>L(R) 



(24) 



The chiral Z' couplings are subjected to constraints from relevant experimental measurements. The 
constraint to B bd has been analyzed 24j, |25(. It is shown 24j that the observed AAIs d and sin 2/3 leads 
to 



y\Re(B L db f\ < 5 x 10" 



(25) 



where y = ( 9 *f, z 



, if only left-handed couplings are considered. Taking y ~ 10 [24|, one has 



B d L b ~ 0.7 x 10- 2 



(26) 



When both left-handed and right handed couplings are included, the constraints are 

y\M(Bi) 2 + [Blf] - 3.8Re(fl&B&)| < 5 x lO" 8 , 
y\lm[(B% b f + (B*) 2 ] - 3.8Im(B&fl&)l < 5 x lO" 8 , 



(27) 
(28) 



which are not as stringent as eq. (|25ll because of the possible cancellation among different terms. 



Using the latest observed AMb s , ref. 25[ shows that 



B 



sb 



10" 



(29) 



To solve the irK puzzle some people have examined new physics in the electrweak penguin sector, 
in particular, the Z' models 3, 27|. Using the experimental data of Br(B — > ttK) in and before 2003, 



one has 



27] 



R c = 1.15 ±0.12, Rn = 0.78 ± 0.10, 



(30) 



which should be roughly equal to one in the SM. It is shown that to explain the deviations a constraint 

'sb D dd 



on B^ b *B dd must be imposed [27 1 



0.01 (solution A) or 0.019 (solution B) 



(31) 
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in the case of £ LL = £, LR . We reanalyze the constraint using the new data of Br(B — > ttK) in this 
year Q . Accoding to the new data, one has [3| , instead of eq. ([501) > 

R c = 1.11 ± 0.08, i?„ = 1.00 ± 0.07. (32) 

With the values, the solution B in ref. [271] is excluded in lcr bounds, but the solution A remains survived. 
So from eas. (|29l3ip . we have 

B& ~ 10. (33) 

The 2T' flavor changing couplings also contribute to B — > 7r°7r°. The constraint from Br(B° — > 
7r°7r°) is: 

0.129 < |0.0887 + 0.0820z + (l.277B^ d - B$ d ) B^\ 2 < 0.178 (34) 

Form eqs. (|26l33p . the above constraints are satisfied. 

Now we obtain the branching ratio of B — > 7r</> in the models for the values of parameters which 
satisfy the constraints discussed above: 

Br{-K<j)) = 1.17 x 10~ 5 |0.0184 + 0.00607i- (B% d + B^) B^\ 2 ~ 10" 7 , (35) 

which is in agreement with the recent data. 

In summary, we have studied the pure penguin process B~ — > 7r~0 using QCD factorization 
approach to calculate hadronic matrix elements to the a s order in some well-known NP models. We 
have shown that the NP contributions in R-parity conserved SUSY models and 2HDMs are not enough 
to saturate the experimental upper bounds for B — > 07r. We have also shown that the flavor changing Z' 
models can make the branching ratios of B — ■+ 4>tt to saturate the bounds under all relevant experimental 
constraints. Therefore, if the data will remain in O(10~ 7 ) in the future it will give a signal of NP effects 
and provide a clue to discriminate well-known NP models. 
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